General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA TECHNICAL 
MEMORANDUM 

NASA TM X-o4943 


ATTITUDE ESTIMATION OF EARTH ORBITING SATELLITES 
BY DECOMPOSED LINEAR RECURSIVE FILTERS 

N75-27036 


Unci ds 
29193 


(NASA-T1-X-64943) ATTITUDE ESTIMATION OP 
EARTH 03BITING SATELLITES BY DECOMPOSED 
LINEAH RECURSIVE PILTP.PS (NASA) 31 p HC 
53.75 CSCL 22A 

G 3/15 


By Shauying R. Kou 
Systems Dynamics Laboratory 


May 1975 


NASA 



George C. Marshall Space Flight Center 
Marshall Space Flight Center , Alabama 


MSFC - Form JI90 (R.v Jun. 1971) 


I NENONT NO. 

NAS/ I'M X-64943 


2. 00' INN Ml NT ACCCSSIC 


I. NECINIENT-S CATALOG NO. 


4 titlc ANO SUBTITLE 

Attitude Estimation of Earth Orbiting Sateilites by Decomposed 


8 REPORT DATE 

May 1975 


Linear Recursive Filters 


6 NINFONMING ORGANIZATION CtCC 


7. AUTMON(B) 

Shauying R. Kou 


M, NC» forming organization NINON r • 


• performing onganization namc and adoress 


10 . WORK UNIT NO. 


George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


11. CONTNACT ON GN ANT NO. 


12 SNONSONING AGCNCV NAMC ANO ADDRESS 


IS. TVNf ON NINON-. & NCNIOO COVENtO 

Technical Memorandum 


National Aeronautics and Space Administration 
Washington, D. C. 20546 


14. SNONSONING AGENCY CODE 


IB. SUNNLEMCNTANV NOTES 

This report is a result of work done by the author while on a post-doctoral NRC fellowship 
at MSFC. 


|«. ABSTNACT 


Attitude estimation of earth orbiting satellites (including Large Space Telescope) 
subjected to t< /ironmental disturbances and noises has been investigated. Modern control 
and estimation theory is used as a tool to design an efficient estimator for attitude esti- 
mation. This new derived estimator is called Decomposed Linear Recursive Filter. It 
is a statistical filter which is optimal in the sense of minimal mean square error. Com- 
parison between this filter and Kalman filter has been discussed. This filter can over 
come the computational difficulty of implementing Kalman filter digitally (integration 
error, roundoff error, etc.). It is shown that simplicity, accuracy and speed are the 
advantages of these filters. Decomposed linear recursive filters for both continuous- 
time systems and discrete- time systems are derived. By using this accurate estimation 
of the attitude of spacecrafts, state variable feedback controller may be designed to 
achieve (or satisfy) high requirements of system performance. 


17 KEY WORDS 

18. DISTRIBUTION STATEMENT 


Unclassified — Unlimited 


✓ ' ^ jr , _s 


/ - ' ' ' si , . . 


19 SECUNITY CLASSIF. (of IhU rapcrtl |2G. SECURITY CLASSIF. (of IhU p«t») 1 2 1 . NO. OF PAGES I 22 PRICE 


Unclassified | Unclassified | 31 I NT1S 

’FC • Fc TO lit! (Rev December 111 !) For tale by Nilional Technical Information S«rvic«. Sprinyfirld, Virtinia 21151 
















ACKNOWLEDGMENT 


The author would like to express his sincere thanks to his 
scientific adviser Dr. Sherman M. Seltzer of the Systems Dynamics 
Laboratory of Marshall Space Flight Center for his discussions and 
help in writing this report. This research was supported by the 
National Research Council. 


TABLE OF CONTENTS 






Page 


INTRODUCTION 1 

SYSTEM DESCRIPTION 2 

STATE ESTIMATION BY KALMAN FILTER 4 

DECOMPOSED LINEAR RECURSIVE FILTER 5 

DI&'l'iCIONS AND CONCLUSION 12 

APPENDIX - DISCRETE-TIME VERSION OF DECOMPOSED 

LINEAR RECURSIVE FILTER IS 

REFERENCES 20 


DEFINITION OF SYMBOLS 


Definition 

transition matrix of a linear system 

transpose of a matrix A 

transition matrix of a composed system 

transition matrix of a linear discrete time system, the 
elements of this matrix are evaluated at time t^ 

transition matrix of a composed discrete-time system 

controlled input matrix 

controlled input matrix for composed system 
controlled input matrix for discrete-time system 
output matrix 
composed output matrix 

output matrix corresponding to state variable x 
output matrix cor; sponding to state variable z 
output matrix for discrete-time system 
output matrix for composed discrete-time system 
matrix related to the dynamics of disturbance z 


the D matrix in the discrete time case 



Symbo 1 


Definition 


K (k) 
x 


K (k) 
x 


K 2 (k) 


related to deterministic disturbances 
related to stochastic disturbances 
the G matrix for composed system 
transformation matrix between variables w and z 
optimal gain matrix of Kalman filter 

optimal gain matrix of Kalman filter for discrete time 
case 

upper part of K^ with n rows 
lower part of K^ with r rows 

optimal gain matrix for a particular initial condition 
upper part of K(t) 



a transformation matrix of Friedland's formula (Eq. (22)) 

a modification of M matrix 

the M matrix for discrete-time case 

power spectral density Function 

covariance matrix 


v 


1 1 

1 1 1 1 ! 

L 

Symbol 

Definition 


P (t) 
x' ' 

covariance matrix with respect to variable x 


p m 

covariance matrix of variable x and z 


XZ' ' 

P z (t) 

covariance matrix with respect to variable z 


P x (t) 

a special l’ x (t) corresponding to a certain initial values 


P 1 

a special case of P(t) 


P 2 

a general form of P(t) 


P k 

P(t) in discrete time case 


P (k) 

X 

P ft) in discrete time case 


P 0 ) 

P (t) in discrete time case 


XZ 

XZ 


P,(k) 

P, it) in discrete time case 





P (k) 
x v ’ 

P (t) in discrete time ca:e 

X 



P(k) 

Pj in discrete time case 


Q 

value of covariance function of s(t) 


\ 

Q in discrete time case 

• 

R 

value of covariance function of v(t) 


R k 

R in discrete time case 


s 

aerodynamic torque and solar pressure torque 

vi 







*k 

T(k) 

T (k) 
x' ’ 

T xz< k > 

V k > 

T(k) 

T (kl 
x' 

u 

U(k) 

U (k) 

x v ’ 

U 2 (k) 


V 



J 1 l 

Definition 

s in the discrete time case 
a posteriori covariance matrix 
the part of T(k) related to x 
a posteriori covariance matrix of x and z 
the part of T(k) related to z 
a special case of T(k) 
the part of T(k) related to x 
control input 

a transformation matrix of Friedland's formula (Eq. (A-16)) 

upper n rows of U(k) 

lower r rows of U(k) 

covariance function of s(t) 

covariance function of v(t) 

sensor noise 

one of the transformation matrix of the formula eq. (22) 
the first two rows of V 
the last three rows of V 



vii 


Symbol 

V 1 

V 2 

v k 

V(k) 

V M 

V z (k) 


w 


X 


X 

X 


X 

*k 

\ 

y 


z 


1111 

Dcf init ion 

modified V form 

another modified V form 

sensor noise in the discrete time case 

transformation matrix for the a posteriori covariance matrix 
upper n rows of V 
lower r rows of V 

gravity gradient torque and earth magnetic torque 

state of a system 

time derivative of the state 

estimated state of the filter 

estimated state corresponding to P x (t) 

state of a discrete time system 

estimated state of 

estimated state of x^ corresponding to f* x (k ) 
output of a system 
y in the discrete time case 
transformed deterministic disturbances 


viii 



Symbol Definition 

9 . estimated val r of z 

z. z in the discrete time case 

k 

z ^ estimated value of 

6 pitch angle of Large Space Telescope 

0 time derivative of 0 

at mechanical angular frequency 

6(t) Dirac delta function 

6^ Kroneckcr delta function 



IX 


TECHNICAL MEMORANDUM X-<i4y43 


ATTITUDE ESTIMATION OF EARTH ORBITING SATELLITES 
BY DECOMPOSED LINEAR RECURSIVE FILTERS 

INTRODUCTION 


The problem of attitude determination (or estimation) of earth 
orbiting satellites is considered in this report. The exactly true 
attitude of those satellites is not available because of the environ- 
mental disturbance torques (e.g. gravity gradient torque, earth magnetic 
torque, aerodynamic torque, etc.) and some undesirable noises (e.g. rate 
gyro noise, CMC tachometer noise, etc.). Hence it becomes necessaiy to 
estimate the attitude of •"te’lites in some optimal way in order to ob- 
tain an estimated (or approximate) attitude and to implement a controller 
by state variable feedback design. It is well-known that Kalman estima- 
tion theory (1) is a useful tool used to estimate the state (or attitude) 
of linear dynamic systems. But a Kalman filter is difficult to be imple- 
mented for higher order systems. Furthermore, if it is implemented 
digitally then the computational accuracy and speed will decrease rapidly 
as the order of the system increases. There are several papers working 
for the reduction of computation burden in the Kalman filter calculations, 
for example, the works have been done by Johnson (2], Simon and Stubheruci 
(3), Samant and Sorenson (4), Friedland (£), etc. The first three of 
these papers are so-called reduced order Kalman filters in which the re- 
quired computations are greatly reduced when only part of the states are 
interesting and being estimated. Friedland's work (5) discussed a com- 
plete state estimation problem. In his paper only the situation of con 

stant bias (or disturbance) is considered. This paper will extend Fried- 
land's work (5) to include the time-varying disturbance case. It will 
be shown that the environmental disturbances and noises on the earth 
orbiting satellites are really either time-varying or random. K'e will use 
the Large Space Telescope (LST), which is an unmanned astronomical obser- 
vatory facility with 3-meter diameter of the primary mirror, as an example. 

It is developed by NASA under the direction of Marshall Space Flight Center 
at Huntsville, Alabama. The description of LST system is shown in section 
2. Special attention is given to model stochastic environmental distur- 
bances and sensor noise as Gaussian stationary white noise processes. For 
the reason of comparison, a Kalman filter for estimating states of the 
LST system is included in section 3. The main result of this paper which 
is called "Decomposed Linear Recursive Filter" (a useful extention of 
Friedland's work (5)) is derived in section 4. In this filter, the covari- 
ance matrix has been transformed into smaller matrices in order to reduce 
the computation burden. Discussions and conclusions arc in section 5. 

A Decomposed Linear Recursive Filter for discrete-time systems has also 
been develope'. Its derivation is slightly different from that of continuous 
time case and is included in the Appendix. References are in section 7. 


SYSTEM INSCRIPTION 


In order to explain dynamics of earth orbiting satellites and some 
physical quantities in detail we use Large Space Telescope as an example. 
The mathematical model of 1ST system described in this section is basically 
according to Schiehlen (6). The LST is modelled as a rigid body having 
reaction wheels as actuators and subjected to gravitational and magnetic 
disturbance torques. These torques arc persistent, deterministic distur- 
bances and can be effectively described by linear differential equations 
(see Johnson ( 7 J ) . In addition, the aerodynamic torque, the solar pres- 
sure torque and the fine guidance sensor noise are also considered which 
cun be modelled as white noise processes. Now for a single-axis analysis 
I.ST system is described, according to (6j, -s 


where 


x 

X 


■ Ax ♦ Bu ♦ Fw ♦ Gs 
2x1 state vector and x 


( 1 ) 

with 0 and § represent the pitch 


angle and the time rate of pitch angle of LST. 


u: a scalar control variable which is related to the driving motor 


torque. 

w: denotes the gravity gradient torque and earth magnetic torque, 

s: denotes the stochastic (or random) aerodynamic and solar pres- 


sure torques. 



The sensor output y corresponding to the pitch motion 8 can be 
represented by 


y ■ 0 ♦ v 


or equivalently 

y « Cx ♦ v (2) 

where C * (1 0] and v is the sensor noise. 


2 


by 


Furthermore, the dynamics of the disturbance torques w can he modelled 


w ■ Hz (3) 

z * Dz (4) 


where z is a 3x1 vector and 


D ■ 


0 0 0 

0 0 2 - 10‘ 3 
0 - 210' 3 0 


H • [1 1 0] 


And the white noise process s(t) of "u. (1) is characterized with 
zero mean and spectral density 1*10"12 arc scc^/scc^. The white noise 
process v(t) of Eq. (2) is characterized by zero mean and spectral density 
8.394- IQ -6 arc s«*c 2 -sec. (These data are obtained from Proise (3)). 


For the following analysis we still need to specify the covariances 
of these white noise processes s(t) and v(t). These covariances can be 
der •••;' from their spectral densities. The covariance function and the 
po « . spectral density function form the Fourier transform relationship 
tor stationary random processes, (see Astrom ( 9] ) . It is assumed that 
white noise processes s(t) and v(t. ) are stationary. Denote the power spec- 
tral density function of s(l) by 

p(u) ■ 1*10"* 2 

for all angular frequency u; then its covariance function 
v 5 (z) V Cov [ s (t ) , s ( t ♦ z ) ] 

■ f*° p(u>) e^ Z du 

•CD 


- 2n-10- 12 6 (z) 

A Q 6 ( z ) (5) 

where 6 ( * ) is the Dirac delta function. Similarly, for the white noise 
v(t), we have 

v y (z) ^ Cov l v (t ) , v (t*z) ] 

• 2n-8.394- 10 -6 6(z) 

- R 6 ( z ) (6) 

Here we introduce notations Q and R for the use of future formula 
tion. Let us rewrite (1) - (4) as a composed system 
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[!] ■ »[:] 


♦ B u ♦ G s 


y » c 


where 


c :i 


C « (C 0] - [1 0 0 0 0] 


Since it is without loss of generality to omit the control term 
Bu in the consideration of estimation problems, let us consider the 
Following simplified system 


[;]-[:]•* 
’ - £ [!] 


where F. ls(t) s (t ) ) * Q 6 (t-t), 
i. [ v (t ) v T (x ) j » R 6 (t-t), 
and s(t ,'(t) are independent. 


STATE ESTIMATION BY KALMAN FILTER 


The Kalman filter is a well-known state estimator for linear 
.. •Hi.vic systems. There are numerous refeiences about the derivation 
.ussions of Kalman filter (e.g. Sorenson 110]). In this section, 
* i show the use of a Kalman filter for state estimation of LST sy- 
When we apply Kalman filter to LST system described by Eqs. (10), 
(11), we have 


Cl ■ • [!] 


* K(t) [y (t) 


■ •- CJ 


where x and 2 are the estimated states of x and z respectively and the 
optimal gain K(t) is defined by 

K(t) « P(t) C 1 R’ 1 (13) 


4 




and P(t) is a SxS covariance matrix which satisfies the following Riccati 
equat ion 

P . A P ♦ P A T ♦ C Q C T • P C T R' 1 C P (14) 

From Eq. (12) we see that for the 5th order system (LST system), 
Kalman filter is also a 5th order dynamic system. The dynamic behavior 
of the covariance matrix P(t) is governed by F.q. (14) which is nonlinear 
matrix differential equation. Furthermore P(t) is a symmetric matrix 
function. 

This is an optimal filter (or estimator) in the sense of minimal 
mean square error between the true state and the estimated stat?. Now 
for the special structures of matrices A, C and G of the LST system 
(see Eq. (9)), the covariance matrix P(t) and the matrix equation (14) 
can be partitioned as follows: 



where Px is the 2x2 covariance matrix of estimation of x. Pz is the 
3x3 covariance matrix of estimation of z. Pxz is the 2x3 cross-covariance 
matrix for estimation of x and z. Then the matrix covariance equation 
(14) takes the form of the following three equations 

P x « A P x ♦ PH P X . T ♦ P x AT ♦ P X2 (FH) T - P x CVIC P x ♦ G Q G T (16) 

f»xz ■ A P i: ♦ FH P z ♦ P xz D t - P x C T R- 1 C P X2 (17) 

P z « D P z ♦ P z D T - P XZ T C T R’ 1 C P xz (18) 

Eqs. (16)- (18) are composed of fifteen nonlinear differential 
equations. These equations are coupled to each other and should be solved 
simultaneously. Furthermore, Eq. (12) is composed of 5 differential 
equations which are coupled to each other too. So to obtain the esti- 
mated state X and 2 we have to solve the simultaneous 15 differential 
equations (Eqs. (16)-ll8)) and simultaneous 5 differential equations (12). 
This is quite computationally time-consuming and has severe numerical 
integration errors. To solve these problems, we derive the Decomposed 
Linear Recursive Filter in the next section. 


DECOMPOSED LINEAR RECURSIVE FILTER 


We would like to point out that the objective of Friedland's 
research [5] is to estimate the state of the following types of systems 

\ x * Ax ♦ Bz ♦ Gs 

I z ■ o (19) 

y *■ Cx ♦ v 
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The second part of this composed system (19), ? » o, represents an un- 
known constant quantity. In this section, the Decomposed Linear Recur- 
sive Filter is derived which is an extension of the above work to include 
the case z ■ D z where D is not a zero matrix, i.e. z is a time-varying 
quantity. Hence this result is applicable to the earth orbiting satel- 
lites including the LST system described by Eq. (10). The basic idea of 
this filter is to use n>3trix transformation of the covariance matrix and 
the particular system structure of earth orbiting satellites to simplify 
the estimation procedure. It turns out that we need to solve only lower 
order matrix equations and estimated states X and 1 can be obtained 
separately. This will greatly increase the computational speed and 
accuracy of estimation of the attitude of satellites. 


At first we investigate a particular solution for the matrix equa- 
tions (16)- (18). Note that Eqs. (I?) and (18) are homogeneous in P xz 
and P z . Hence, if P xz (o) ■ o and P z (o) * o then P xz (t) = o for all 

t > o, then Eq. (16) becomes 

P x « A P x ♦ P x A 1 ♦ G Q gT - P x C T R" 1 C P x (20) 


the notation P x is used to poin' out tue ^orms of Eqs. (16) and (20) arc 
different. So for this particular selectman of initial values of covari 
ance matrix, let us denote the solution of Riccati equation (14; by 


Pj(t) 


P (t) 0 

x v ' 

0 0 


( 21 ) 


The relation between solutions of Riccatt equations of the form 
(14) corresponding to different initial coudit'.ons can be expressed 
according to a result of Friedland [11]. 


Theorem ([11]): If P (t) is a solution to Riccati equation (14), any 
other solution P,(tj of Eq. (14) correspond ng to a different initial 
condition can be"expressed as follows: 


r 2 (t) = Pj(t) ♦ VMV 


( 22 ) 


-T d-1 


where V is a 5x3 matrix satisfying V - (A - Pj C‘ R 1 C) V (23) 

and M is a 3x3 matrix satisfying M = -M R _ * C V M (24) 


and initial condition of P 2 (t) is 


P 2 (o) 


P (o) 
x v ’ 


P (o) 
xz v ’ 


P xz V°) 


(25) 


where P x (o) is the nonzero part of the initial condition for Pj(t). 

In this theorem, we see that M is a symmetric matrix and V can be 


6 



J I I ...i 


partitioned as follows 


V - 



(26) 


where is 2x3 matrix and is 3x3 matrix. Then by using the definition 

of A and Pj, Eqs. (23), (24) can be expressed by 

V x - (A - P x C T R* 1 C) V x ♦ FH V z (27) 

V z - D V z (28) 

M - -M V x T C T R' 1 C V x H (29) 


Now let P,(t) be partitioned according to (15), then 
P ,<‘> ■ F x<*> ♦ V x M < 

P XI (,) ' v x M v x T (30) 

P ,<‘> - V z M V, T 

So the initial conditions V (o), V (o), M(o) of Eqs. (27)-(29) must 

satisfy 

P x (o) « P x (o) ♦ V x (o) M Vj (o) 

P xz (0) " V x (o) M(o) V 7 (o) 

P z (o) - V z (o) M(o) V z T (o) 

If P (o) is also assumed to be zero, then we have 
xz v 

V x (o) * o 

V z (o: - I (31) 

Mv.j) = P z (o) 


We want to show that the inverse of the square matrix V exists. 
From Eqs. (28) and (31) we have 

V z (t) * e Dt V z (o) - e Dt 

which is nonsingular (see Brockett (12]). Hence its inverse V, 1 exists. 
Now for the purpose of simplification of equations used for estimation of 
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: 


states of LST system, we introduce new variables Vj and V 7 

-1 

I V V 

-1 


Vj ■ v v 2 _1 


V 2 ■ V x V‘ 


V x v x 


(32) 


(33) 


Here we have V. is also a 5x3 matrix consisting oi two parts-- 
* 

the lower part is a constant identity matrix. Again we introduce new 
variable M, 


’I 

M 1 ‘ V z M V z 


-1 


(34) 


Now we try to express P z ( t ) in terms of these new matrices Vj, V 7 
and Mj. From Eq. (22), we have 


P 2 (t) * P i * V i V z M < V l V 

' P . * V 1 M 1 V i T 


(35) 


And from Eq. (30) we have 

p x (t) ■ P„(t> <■ V x »,-* H » t T ,V 2 T )-‘ V x T 
■ V * 1 * V x V' \ M V x T < V x' 1)T V x T 


P x (t) * V, Mj V 2 ' 


(36) 


P (t) * V V _1 V M V T 
XZ V ' X z z z 


V 2 Mj 


(37) 


P 2 (t) - V z M V z 


M, 


(38) 


So far we know that to obtain P . P v , and P we need only to have 

X a z z 

V.,, M, and P . P can be obtained from Eq. (20). Now we will show that 
2* 1 xx 

V 7 and Mj satisfies the following equations: 


V 2 =« (A - P x C T R" 1 C) V 2 ♦ FH - V 2 D 


(39) 
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(40) 


1 


Hj ■ D Hj - Mj V 2 T C 1 R' 1 C V 2 Mj ♦ M J I) T 


-1 


Proof: Since m V x^z • 


• • -1 • -1 

V, ■ V V 1 ♦ V (V ) 
2 x z x v z ’ 


V V _1 - V V ** V V _1 
xz x z z z 


- [ (A - P C T R" 1 C) V, ♦ FH VI V/ 1 - V V “* D V V, 

A A •• v A t V « 


-1 


-1 


-1 


(A - P C T R' 1 C) V, ♦ FH - V, D 
v x '2 2 


Similarly M. ■ V M V , so 
l z z 


• • 


M-V MV T *V M V T ♦ V MV T 
1 z z z z z z 


D V MV T + V (-M V T C T R" 1 C V M) V T ♦ V M (D V ) T 
zzz'x x ’ z z ' z 1 


DM 1 - Mj V 2 r C T R 1 C V 2 Mj ♦ Mj D 1 


Remark: In the above analysis, we see that to find covariance matrix l’-,(t) 

(which is governed by a nonlinear coupled matrix differential equation (14) 
and should be solved simultaneously as shown in section 3), we need only 
find out V 2 (t), Mj(t) and P x (t). Now the important thing is that the 

equations ((20), (39) and (40)) of V,, Mj and P^ are not coupled in the 
following sense. We can solve P x by using Eq. (20) without using any in- 
formation about V 2 and . Next we obtain V 9 from Eq. (39) and information 
P x> Finally using V 2> F x and Eq. (40) we obtain Mj . As a matter of fact 

these matrices can be solved separately. It is obvious that by using the 
method of this analysis to find covariance matrix, we can save a lot of 
computation time and increase greatly the numerical integration accuracy. 

Next if we substitute the results Eqs. (36) - (38 ) into Eqs. (12) 
and (13) we have 

i * a « ♦ fh 2 ♦ (P x ♦ v 2 M ! V 2 T ) cT r1 (y - C *) (41) 

i = D 2 ♦ Mj V 2 T C T R _1 (y - C 2 ) (42) 

Now we want to show that how & and 2 can be obtained separately 
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Theorem: The solution X of Eq. (41) can be expressed by 

X • x ♦ V 2 2 

where x is the solution of the following system 
X ■ A x ♦ F x C T R ' 1 (y - C x) 
and V.,, 2 satisfies Eqs. (39) and (42). 


(43) 


(44) 


Remark: x(t) can be obtained from Eq. (44) whenever we have l’ x and y(t). 

It is in fact the estimate of state x(t) when there is no deterministic 
external disturbance torques, i.e. z » o. 

Proof: To prove this theorem, we have to check that the expression (43) 

satisfies Eq. (41) or equivalently the following formula should be equal 
to zero 

X - A X - FH 2 - (P x ♦ V 2 Mj V/) C T R’ 1 (y - C X) - 0 (45) 

In the following we will substitute all x by x ♦ V-, 2 to the left 
hand side of Eq. (45). First we give an expression for x, 

5 * * * Ht V 2 2 
* x ♦ V., 2 ♦ V, £ 

- A x ♦ P x C T R' 1 (y - C x) M(A - P x C T R" 1 C) V 2 ♦ FH - V , D) 

♦ (D 2 ♦ M V, T C 1 R' 1 (y - C (x ♦ V, 2))) 

Now the left hand side of F.q. (45) becomes 

A x ♦ P x C T R" 1 (y - C x) ♦ t(A - P x C 1 R" 1 C) V 2 ♦ FH - V 2 D) 2 

♦ V 2 (D 2 ♦ Mj V/ C T R _1 (y - C (x ♦ V 2 2))] - A (x ♦ V 2 2) - FH 2 

' (R x * V 2 M 1 V 2 T) cT Rl (y " C (* ♦ v 2 2)) 


(- P x C T R" 1 C V, - V 2 D) 2 ♦ V ; [D 2 ♦ Mj V 2 ‘ C* R' ‘ (y-C(x*V., 2))) 


T _T 1 


- V 2 M l V 2 T C T R 1 (y - C(x ♦ V 2 2)) ♦ P x C R _ * C V 2 2 


- T -1 

n f * 
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This completes the proof. 

Remark: Based on this theorem and the above analysis, we will summarize 

the procedures in finding the estimated states 8 and 2 as follows: 

(a) From Eq. (20) we have F . 

(b) Using P x and Eqs. (39) and (40) we get and M . 

(c) Again using P x and the system's output y(t) we get x(t) from 

Eq (44). Note that in this step we don't use any information 
about 8 oi 2. 

(d) Now wc use the results of (a) - (C) to find 2 without using 
any information about 8. That is frori Eqs. (42) and (43). 

2 - D 2 ♦ Mj V, T C T R" 1 (y - C (x ♦ V 2 2)] (46) 

(e) At last from Eq. (43) we I *ve 

8 - x ♦ V 2 2 

In this procedure we see that the estimation of states of LST 
system has been decomposedinto two parts; 2 and 8. Estimation of z can 
be obtained by using data x instead of 8. Furthermore x and 2 are ob- 
tained from Eqs. (44) and (46) which are 2nd-order and 3rd-order systems 
respectively. On the contrary, if the Kalman filter is used, we have 
tc build up a fifth order dynamic system (filter) to estimate the states 
of the LST system. So by using the analysis and techniques of this sec- 
tion, filters can be implemented to obtain estimated states 8 and 2 by 
lower order systems. In addition, the simplicity in structure of this 
Decomposed Linear Recursive Filter gives us considerable amount of re- 
duction in computation burden and storage requirements. 
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DISCUSSIONS AND CONCLUSION 


(a) Discussions 

1. High order estimators usually are too complex for practical 
mechanization and are very parameter sensitive. Now if we use Decomposed 
Linear Recursive Filter to estimate the attitude of LST system, we need 
only to build up 2nd-order and 3rd-order dynamic systems (instead of a 
Sth order Kalman filter) as an estimator. 

2. If the estimators or filters arc implemented digitally, Che 
computations required by a high-order filtering algorithm may become 
excessive (in memory capacity, accuracy and computational speed). 

J. The characteristic of this ;uw filter is that the low order 
covariance matrix Eq. (20) is independent of matrices V and M of Eq. |22). 
Similarly, the 2nd-order system X can be obtained from Eq. (44) independent 
of the value of 2. This property is a kind of decomposition of a set of 
equations and greatly simplifies the estimator design. 


(b) Conclusion 

In this paper, a Decomposed Linear Recursive Filter has been 
derived for the estimation of attitude of earth orbiting satellites in 
wnich the disturbances and >iscs are either time-varying or random. Sim- 
plicity, accuracy and speed are the advantages of this filter. This filter 
ran be easily implemented by the on-board digital computer and give an 
accurate estimation of attitude of satellites. Furthermore; this is the 
necessary information for state variable feedback controller design. 

The analysis of section 4 can be generalized to linear time-varying 
systems such as 

| x ■ A(t) x ♦ B(t) z ♦ s 

I z « D(t) z 

y « Cj x ♦ C, z ♦ v 


where A, B, 1), Cj and C 7 can be matrices with time-varying elements, and 

y may also include the disturbance z explicitly. The proof of this fur- 
ther extention is exactly the same as that of section 4. 
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APPENDIX 

Discrete-Time Version of Decomposed Linear Recursive Filter 


Most of the practical implementations of filters are by digital 
means. Here we derive the analog counterpart of Decomposed Linear Recur- 
sive Filter for discrete-time systems. 

Consider 

x k ' \-i Vi * Vi 'k-i * Vi (A - n 

'k ■ D k-i 'k-i '*-*> 


where 


x^ is the value of n-dimension state at time t ^ . 

z, is the value of r-dimension deterministic environmental dix- 
ie 

turbances at time t. , 

k 

\ 1* B k 1’ 1 are time * var y in K matrices with dimensions n x n, 

n x r, r x r respectively and assume that D^ is nonsingular for 
all time t 


k* 

s^ j is n-dimension white noise sequences with 
E (*k *j T l " Q k 6 kj 

where 5. . is the Kronecker delta function. 

kj 

and the measurement equation is expressed by 


c k x k * v k 


where 


y^ is the value of m-dimension output at time t^, 

C. is a m x n time varying matrix, 

v^ is m-dimension white noise sequences with 

E ' v k v j t i ■ *k s kj 


(A-3) 
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t.ombinc (A- 1 ) and (A- 2) , we have 


• Vi 


and 


Sk 


k-1 


k-1 


♦ v. 


♦ G s 


k-1 


C A - 4 ) 


(A-S) 


where 


*k-l 


Vi 

0 


k-1 


k-1 


C -k 


(C k 0) 


Kalman filter formulation for system (A-4), (A-5) is 


*k-l 


k-1 


k-1 


h l y k " ^k - 


^-1 


‘k-1 


k-1 


(A-6) 


where i. and 2. are the estimated values of x. and 2 . , 
k k k k 

K ■ Pb C. 1 (C. P. C. T ♦ Rj’ 1 
k k — k l -k k -k k J 


(A-7) 


Py is the a priori covariance matrix satisfying the following 
difference equation 


Vi ' 4k " - “k V p k V * c cT 


( A- 8 ) 


Let 


P x 00 

p T (k) 
xz v ' 


P (k) 
xz' ' 

P z (k) 


(A-9) 


where 

P (k) is n x n covariance matrix of estimate of x at time t. , 

X ^ 

P _ ( k ) is r x r covariance matrix of estimate of z at time t^, 

P x ,(k) is n x r covariance matrix of estimates of x and z. 

According to the submatriccs of Eq. (A-9) and using ( A - 1 ) . (A-2) and (A-3), 
we have Eq. ( A- 8 ) being expressed by the following three equations 
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I 




P.0..1) • c \ Cl - < V * »k K c k' I p ,< k > V ♦ p „ (k) *k T i 


• \ I \/< k > \ T ♦ p l' k > »k T I • «k.l 


(A-10) 


p ,i (k * n • i \ c • s* c k> * B k k k‘ c k i p « t (k » D k T * , k p . ,k) V 
P ,< k *'> ■ - D k C k P xx< k > D k T * D k v k > °k T 


where 


K k * 


X z 

, is a (n ♦ r) x m matrix and are n x m 


and r x m respectively. 

If P (o) ■ o, P (o) ■ o, then we have the solution of (A-8) denoted by 
Z X z 


P (k) 0 

x* ' 


(A- 13) 


0 0 
for .ill k >o and P x (k) satisfies 

p x (k4l > • I \ - <\ k k ♦ 8 k *k‘> C k 1 F x (k > V * ’k.l (A -‘ 4) 

By Kalman filtering theory (see Sorenson (10)), a priori covariance 
matrix P can be related with a posteriori covariance matrix T by 

X * 


P x (k-1) . A k T x (k) A k * ♦ Q k4l 


(A- 15) 


Note: By definition, a priori covariance matrix l‘ x (k ) is defined by the 

covariance matrix of the estimate of x(t k ) given the observations 
{ y l p y 2* *••• y k but not y k‘ A P ostericTi covariance matrix T x (k) 
is the covariance matrix of estimating x ( t k ) given the observation 
| yj, y,, ...» y k ^ , i.e. including the current output y k> 

Theorem (Friedland (11)) If P (k ) is a solution of equation (A-8), then 
any other solution P(k) of Eq. (A-8) corresponding to different initial 
Conditions can be expressed by 


P(k) - P(k) ♦ U(k) M(k) U'(k) 


(A- 16) 


where 


U(k) is a (n ♦ r) x r matrix and M(k) is a r x r matrix satisfying 
U(k*l) - \ ( I - POO Cj (^ P(k) cj* R k ) _1 ^ ] U(k) ( A- 1 7 ) 
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H ( k ♦ 1 ) - H(k)-M(k) U T (k) cJiC^ P(k) cj ♦ R k ♦ ^ U(k) M(k) U T (k) 1 ] ' ' 


• U(k) M (k ) (A- 18) 

In terms of this transformation, (A- 16) can be rewritten as 

P(kM) • F(k*l) ♦ U (k ♦ 1 ) M(k^l) U T (k*l) (A- 19) 

- P ( k ♦ 1 ) ♦ ^ [ I • K(k) C fc ) U(k) M(k*l) U T (k) (I - K(k) C R J 1 aJ 

(A- 20) 

where 

K(k) - P(k) C, T [ P(k) cj ♦ R k I* 1 (A-21 ) 


Similarly, we have 

P(k*l) - ^ T(k) A^ ♦ G g kM G T (A-22) 

P(k*l) - ^ T(k) A Jt T ♦ G Q kM G T (A-23) 

Hence (A-20) can be expressed as 

T(k) - T(k) ♦ V (k) M(k*l) V T (k) (A-24) 

where 

V (k ) ■ [ I - K(k) ] U(k) (A-2S) 

and by (A- 17) we have 

U(k*l) ■ ^ V (k) (A- 26) 


Using (A- 13), Equation (A-21) turns out to b 



’V k > C k T 1 C k " P x (k) C k T * R k i' 1 

A 

K (k) 

X 

K(k) ■ 

0 

m 

. 0 . 


Let U(k) and V (k ) be also partitioned as follows 



U (k)‘ 


V (k) 

U(k) - 

A 

. V(k) - 



u z (k) 


M 

> 

1 


where is n x r, U, is r x r, V is n x r and V. is r x r. 
Then, using (A-27), V(k) and U ( k ♦ 1 ) of (A-25), f A- 2ft ) become 

j V x (k) * ( I - K x (k) C k ] U x (k) 

I v z (k) - u 2 (k) 


C A - 28 ) 


(A-29) 
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| V k * l > ■ \ V k) ♦ B k V k > 

I U 2 (k>l) - D k V z (k) 

Finally, F.q. (A- 19) become* 

P x (k) - P x (k) ♦ U x (k) M(k) u; T (k) 

P x: (k) * U x (k) M(k) U z T(k) 

P z (k) - U z (k) H(k) U z T (k) 

and Eq. (A- 24) becomes 

T x (k) - T x (k) ♦ V x (k) M(k*l) V x T (k) 

T xz (k) - V x (k) H(kM) V z T (k) 

T z (k) - V z (k) M(k*l) " z T (k) 

It can be shown that " of (A-7) can be expressed by 

*k • T(k) cj R / 1 


•c x 00 

K z (k) 


T x C k T R k _1 
T xr T C k T "k 


-1 


Substituting T x 


and T 


xz 


according 


o (A-32) , we have 


K x lk) * (T x (k) ^ V x M(kM) V x T) C k T R k’* 
K z (k) - V z (k) M(k*l) V x T (k) C k T R k ** 


(A - 30) 


(A- 3 1 ) 


(A-32) 


(A-33) 

(A-34) 


(A-3S) 

(A-36) 


Now from (A-6), 

Vl "Sc *k * B k h * V* 1 ’ 1^.1 - C k.! \ 8 k - c k.| B k *k* (A 


and 


2. , ■ K K ♦ K (k*l) l y t , • C. . A. i • C. . B, i. ] (A- 38) 

k-*l k k z v ’ 1 7 k*l k>l k k k*l k k 1 


Theorem: The solution of Eq. (A-.’’) can be expressed as 

2. ■ x. ♦ V (k) V/ l (k) 2. 
k k x z ' k 


(A -39) 


-37) 
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where if, the solution of the following equation 


and 


Vl “ '*k X k * K x (kM) * y k*l " C k*\ \ X k * 


M k > ■ ? .< k > C k T I C k V k > C k T * R k I'* 


(A-40) 


and V , V and 2. satisfying (A-29) and (A- 38 ) respectively. 

X Z K 

Proof: Wc first show that V z * exists. Since 

V z (kM) - D k V 2 (k) 

V 2 (o) - I 


so 


v 2 (k-l) . D k D k-l 


is not singular rince are assumed to be nonsingular as in ( A - 2 ) . So 
V, 1 exists. 

Next we wilt prove that both sides of F.q. (A-37) are equal by substituting 
the expression (A-3. ) into Eq. (A-37). 

Left hand side of F.q. (A-37) becomes 

X. , * x. , ♦ V (k*l) V - 1 (k*l) S.. 
k*l k*l x z k*l 


-1 


\ «k * V k “> &k.l - c k.l A k *k> * V , (k *‘> V x (k4,) 


[D k t k . *,».» (y k „ - C k>1 A k S k - C k>1 B k t k ) 1 

■ A k ; k * V k4 » l^k.l - C k.l A k \ ] * V k4l > V 1(k4l) 

[D k S k * K,(k*l) (y k ., - c kM A k r k - c kM (A k k k v ^ 1 . 
B k ) 2 k ] (A-41) 


Right hand side of Eq. (A-37) becomes 
\ »k 4 B k 2 k 4 k x (k4,) 1^1 - 

- \ J k 4 \ v , (k > V‘ (k > 5 k 4 


C k*l A k x k " C k*l B k \ 1 

B k \ + K x (kM) l) W C k*l A k VSflVk 1 
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\ v \ v « v «'‘ 5 k * \ v K . (k *" 'mi-Mi www; 1 - 


v ! k > 


(A-42 ) 


By (A-33) and (A-32) wc have 

K x (kM) “ [ V k<1) * V** 1 ’ M(k * 2 > V x T ( k ^^ ) C Jl Vl 

■ V k * n c Ji Mi ♦ v ,< k *‘> M < k * 2 > v x T < k *') c Ji Ml 


K x (k *l) ♦ V x (kM) V z * I (k*l) K z (kM) 


(A-43) 


Comparison of expressions (A-41) and (A-42), and using the relation (A-43), 
it turns out we need only to prove that 

\ v x (k) V 1(k) ♦ b k - V k *n c vi ( \ v x (k) v 2 " 1(k) 4 B k ) * v x (kM) 


V 1(k M) D k - V x (k*l) V z - 1 (k.l) K z (kM) C kM (A k V x (k) V 2 - 1 (k) ♦ B k 1 

( A- 44 ) 


Proof: left-hand side of (A-44) 

■ S V k > V'< k > * B k - V k *’> Ml (A k V x< k > V'< k > * “k> 

- V x (k.l) V s -‘(k.l) k,(k*l) C k „ (* k v x V,-* . B k ) 

- ( 1 - K x (k*l) C k+1 ) (A k V x V z _1 ♦ B k ) - V x (k>l) V z * , (kM) K z (k*l) 

Ml lA k V x V x'' * B k> 

. ( I - K x (kM) C k .,) V x . B k V,) V 1 - V x (k.l) »,- l (k*l) 

M k *» C k4 , (* k V x V I ' 1 . B k ) 

- v x (k-i) v I - , (k) - v x v t -* k 2 c k „ v x » t -' • B k ) 

- V x (k.l) V 2 -'(k.l) 0 k - V x (k.l) V 2 -‘ k 2 C kM (A, V x (k) V 2 -'(k) • B k l 
= right hand side of Eq. (A-44). 

This completes the proof. 
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